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ABSTRACT 

There is increasing evidence of a connection between AGN activity and galaxy evolution. 
To obtain further insight into this potentially important evolutionary phase, we analyse the 
properties of quasar host galaxies. In this paper, we present a population synthesis model- 
ing technique for off-axis spectra, the results of which constrain host colour and the stellar 
ages of luminous quasars (Mv{nuc) < —23). Our technique is similar to well established 
quiescent-galaxy models, modified to accommodate scattered nuclear light (a combination 
of atmospheric, instrumental and host galaxy scattered light) observed off axis. In our model, 
subtraction of residual scattered quasar light is performed, while simultaneously modeling the 
constituent stellar populations of the host galaxy. The reliability of this technique is tested via 
a Monte-Carlo routine in which the correspondence between synthetic spectra with known pa- 
rameters and the model output is determined. Application of this model to a preliminary sam- 
ple of 10 objects is presented and compared to previous studies. Spectroscopic data was ob- 
tained via long-slit and integral-field unit observations on the Keck and WIYN telescopes. We 
confirm that elliptical quasar hosts are distinguishable (bluer) from inactive ellipticals in rest 
frame B-V colour. Additionally, we note a trend for radio luminous {L^ghz ^ 10^*^ erg s^^) 
quasars to be located in redder host galaxies in comparison to their less luminous radio coun- 
terparts. While the host colour and age of our radio luminous sample is in close proximity 
to the green valley, our radio faint sample is consistent with quiescent star-forming galaxies. 
However, further observations are needed to confirm these results. Finally, we discuss future 
applications for our technique on a larger sample of objects being obtained via SALT and 
WIYN telescope observing campaigns. 
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;> 1 INTRODUCTION 
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It is well established that active galactic nucleus (AGN) activity is 
due to accretion onto super massive black holes (SMBHs). Addi- 
tionally, SMBHs are found to be common, existing in most, if not 
all, massive galaxies. Observationally, the SMBH mass is found to 
be tightly corr elated to the bulge ve locity dispersion raised to the 
fourth power ( lGebhardtet"ai]|200(]|) . However, the bulge extends 
well outside the gravitational influence of the SMBH. This correla- 
tion has led to t he hypothesis that the bulge and SMBH co-evolve 
( lKormendvll2000,l Additionally, models of galaxy evolution must 
include regulation mechanisms to adjust and quench star formation 
in massive galaxies after early times to prevent the over production 
of bright galaxies ( Benson et al. 2003). 

Evidence pointing toward the co-evolution of SMBHs / bulges 
and the need for regulation mechanisms has led to models which 
utilize AGN feedback to quench star formation. An example of a 
m odel of relevance to ou r studies is the hydrodynamic simulations 
of iHopkins et al. I ( l2006h . In this model, gas rich mergers induce 
gas inflow triggering both star formation and quasar activity. Some 



small percentage of AGN energy output is converted to thermal 
energy, assisting in the stoppage of further star formation and ac- 
cretion. Simulations, which feature self-regulated SMBH growth, 
have been successful in reproducing variou s observables, such a s 
the bimodal colour-magnitude distribution dCattaneo et alj|2006h . 
However these models must operate on a very large range of scales 
(fipc - Mpc) and must incorporate many physical processes. There- 
fore it is desirable to further constrain these models observationally. 
The goal of this work is to develop an off-axis method to provide 
stellar age and host colour constraints during the quasar phase. 

Unfortunately, studies of the underlying stellar populations 
are hampered by the overwhelming emission from the cen- 
tral quasar. Consequentially, progress has been slow and, in 
some cases controversial, in the struggle to place observational 
const raints on quasar h o st galaxies. One g r oup of collabora- 
tors ('McL ureetaQligggl ; iHughes et al.l l200d : iNolan et al.l 1200 ll : 



[Punlop et air 2003 ) believes that local quasar (Mv{nuc) < 
—23.5) hosts are predominantly normal massive elliptical galax- 
ies. However, Hubble Space Telescope (HST) morphology 
( iHooper. Impev & Folt j 1 19971 ; iBennert et al] l2008h . multicolour 
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Object name redshift Mv(nuc) My (host) 



logiLscHz) Host 
ergs^^ Morphology 



3C273 0.1583 -26.7 -23.2 44.1 E 

4C 31.63 0.2950 -25.1 -25.1 43.3 E 

PKS 1302-102 0.2784 -25.9 -22.9 43.0 E 

PKS 0736+017 0.191 -23.2 -22.6 43.0 E 

PKS 2135-147 0.2003 -24.7 -22.4 42.9 E 

PKS 2349-014 0.1740 -24.5 -23.2 42.5 E 

PG 1309+355 0.1840 -24.4 -22.8 41.3 S 

PHL909 0.171 -24.1 -22.2 40.0 E 

PG 0052+251 0.1550 -24.1 -22.5 39.4 S 

PG 1444+407 0.2673 -25.3 -22.7 39.2 S 

Table 1. Test sample properties. My (nu c) is the V-band abso lute magnitude of th e quasar. My (host) is the V-band absolute magnitude of the host galaxies. 
Both of these quantities are derived from lBahcall et alj h997h when possible, see IWolf&Sheini3<2008ll for d etails. Radio luminosity i s estimated utilizing 
the NASA/IPAC Extragalactic Database and assuming spectral index of -0.5. Host morphology is obtained via lHamilton. Casertano &"Tum shek I2002J. "E' 
denotes elliptical and 'S' denotes spiral. 



im aging (jjahnke. Kuhlbrodt & Wisotzkill20 04l). and spectral stud- 
ies jCanalizo & Stocktonl200ll : lMiller &S heinis 2003|) have found 



evidence inconsistent with a population of normal relaxed ellipti- 
cals. In fact, an on-axis spectroscopic study probing the inner re- 
gion of nearby quasar ho sts found young Sc- hke stellar populations 
in half of their sample tLetawe et al.ll2007h . The discrepancies in 
these studies may be due to relative ly sma ll sample sizes, system- 
atics, and the different radii probed (Lacy '2006). Additional study 
is required to develop a clear understanding of stellar properties of 
nearby luminous quasar hosts. 

In IWolf& SheinisI j2008l) we presented velocity dispersion 
measurements of a sample of nearby luminous quasar host galax- 
ies. Future work will examine the relation between SMBH mass 
and velocity dispersion (Sheinis et al. 2010 in preparation). In this 
paper, we develop a method to measure the stellar age and colour 
of nearby quasar host galaxies at off-axis radii of 9 to 15 kpc. We 
observe off-axis spectra to minimize the observed quasar emission, 
thereby maximizing signal-to-noise. However even in off-axis ob- 
servations, described in §2, scattered light from the central source 
can still contribute significantly to the observed light. Our tech- 
nique of removing this scattered light is similar to previous stud- 
ies (Boroson, Persson & Oke 1985; Miller & Sheinis 2003) and de- 
scribed in detail in §3.1. The modelling of the host spectra (§3.1) 
and the error es timation technique (§3.2) are based on qui escent 
galaxy models jTremontil |2003| ; ICid Femandes et ak II200I here- 
after CF05). Application of our model to our preliminary sample 
is found in §4. A comparison to recent spectroscopic and photo- 
metric studies is found in §5. Summary and conclusions are pre- 
sented in §6. Appendix A contains comments on individual objects. 
Throughout this paper a cosmology of Ho ~ 70 km Mpc~^ , 

~ 0.3, and Qa = 0.7 is adopted. 



2 OBSERVATIONS 
2.1 Sample 

Our ongoing spectroscopic observing campaign has resulted in 
a current sample of 28 nearby luminous quasars. This sampl e 



discussed below. Additionall y, this subset of ten obj ects has known 
stellar velocity dispersions ( IWolf & SheinisI |2008|) . eliminating a 
potential free parameter of the model. This relatively small sam- 
ple of quasars is not representative of the local quasar population, 
which consists of ~10% radio loud quasars. Our sample contains 
six radio loud and four radio quiet objec ts. We define an object a s 
radio loud by the criteria established bv lKellermann et al.l ( Il994l) . 
L5 GHz > 10^^ W Hz-^, which is roughly lO^^'^erg s'^ for our 
adopted cosmolog)0. All quasars are local z < 0.3 and luminous 
Mv < —23. The properties of our test sample are summarized in 
Table 1. 



2.2 Spectra 

Data ac quisition and red uction is des cribed in detail in SheinisI 
( |2002|) . iMiller & Sheiiiii (^003), and IWolf & Sheinii ( |2008|) . In 
summary, seven of the ten object s were observed via the Low 
Resolution Imaging Spectrograph ( IOkeetal.|[T994l) on the Keck 
10-m telescope. The remaining three objects we re observed us- 
ing a 82 fiber integral field unit (IFU), SparsePak teershadv et al.l 
2005) which feeds the Bench Spectrograph on the 3.5-m WIYN 
telescope. Short on-axis observations, typically 1-2 minutes for 
Keck and 30 minutes for WIYN, and longer off-axis (2"-4.5" from 
the nucleus) observations, typically 0.5 - 4 hours, were obtained 
for each object. The approximate off-axis slit and fiber positions 
with respect to archival imaging is displayed in Fig. 1. The prop- 
erties of the off-axis pointings are summarized in Table 2. Keck 
objects have a rest wavelength range of ~3500-7000A at a spec- 
tral resolution of A A ~ 11 A. WIYN objects have a rest wave- 
length range of ~3500-6000A at a spectral resolution of AA ~ 
5A. All spectra were corrected fo r Gala ctic extinction using the 
law of Cardelli, Clavton & Mathi^ (| 19891) and the Ay values from 
ISchlegel. Finkbeiner & Davis. 419981) as listed in the NASA/IPAC 
Extragalactic Databas43 (NED). 



consists of ob ects previously imaged by ipancaii ei ai.i iiyy/v . , r— -1 r-rr-i , . „ , 

.,, r r 1. /• |_, , ^^n~nAr.V.I , IKellermann et al.| U 994ft observationally determined the luminosity 

With the addition of a tew obiects from IDunlop et all 120031) and »u 1. u • 1 r rr r,r,, _i , , -1 -, m 



Bahcall et al 



jl997h 



iOuvon. Sanders & Stockt on (2006). We demonstrate the capabili 
ties of our spectral synthesis model on a subset of 10 objects. These 
objects have spectral data that display sufficient signal-to-noise and 
stellar continuum to constrain the stellar properties via the method 



threshold assuming a cosmology of Ho = 70kms~^ Mpc~^ , go = 1/2. 
We solve for the observed flux at z = 0.2 and then convert to our adopted 
cosmology. The specific luminosity is multiplied by the observed frequency 
to obtain the quoted 5 GHz luminosity. 
^ http://nedwww.ipac.caltech.edu/ 
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Figure 1. Sample of observed quasar hosts with approximate off-axis slit and fiber positions indicated in red. The slit lengths are much lai'ger than indicated, 
extending approximately 7' projected on the sky. For WIYN objects, only the off-axis fibers used for the host galaxy spectrum are shown. Objects are listed 
from left to right, top to b ottom by most radio bright to least. All images are from HS T (Bahcall et al. 1997) (23"x23") with the exception of 4C 31.63 
(10"xlO") jGuvon et al.]|2006) and PKS 0736+017 (30"x3q ") fautchings. Johnson & Pvke 1988) . Where indicated, the 10 kpc scale is computed for 
Co = 10 and Hq = 100 km s"! Mpc'^ jBahcall et aljl997h . For our adopted cosmology this scale corresponds to ~ 15.7 kpc. 



Pointing 


Robs (kpc) 


S/nA-1 
(5500-5700 A) 


Quasar Scattering 
(3600-5900 A) 


Telescope 


3C273 4N 


11.79 


18.6 


68.8% 


Keck 


4C 31.63 2N 


8.74 


11.5 


50.5% 


Keck 


4C 31.63 2.5E 


10.92 


9.4 


34.6% 


Keck 


4C 31.63 3S 


13.11 


10.1 


38.7%c 


Keck 


PKS 1302-102 2.3N 


12.79 


9.3 


68.4%o 


Keck 


PKS 0736+017 4.5NW 


14.18 


8.3 


49.5% 


WIYN 


PKS 2135-147 3W 


12.47 


12.3 


70.8% 


Keck 


PKS 2349-014 4N 


14.12 


14.8 


21.8% 


Keck 


PKS 2349-014 3S 


9.47 


18.8 


25.2% 


Keck 


PG 1309+355 4.5SW 


13.76 


13.4 


41.3% 


WIYN 


PHL 909 4.5N 


12.97 


11.7 


77.8% 


WIYN 


PG 0052+251 3S 


9.65 


23.0 


64.6% 


Keck 


PG 1444+407 3S 


15.07 


8.8 


69.6%, 


Keck 



Table 2. Off-axis spe ctral properties. Pointin g designation indicates object and approximate off-axis position in arcseconds. ifobs is the observed off-axis 
radius, as reported in IWolf & SheinisI )2008h . Additionally, off-axis signal-to-noise, fraction of the observed off-axis spectrum which consists of scattered 
light, and telescope are shown. The scattered quasar light percentage is determined by the model output. 



3 SPECTRAL SYNTHESIS 



3.1 Model 

Our spectral s ynthesis code is a modified version of the model first 
discussed in .Miller & SheinisI ^2003), with the underlying stellar 
spectrum now modeled as a weighted summation of 15 instanta- 
neous starbursts of various post burst ages. We model the off-axis 
spectrum, which consists of a stellar component and a scattered 
quasar component. Quasar light is scattered into the off-axis line of 
sight due to atmospheric seeing effects, as well as a small contribu- 
tion from the intrinsic optical point-spread function and astrophys- 
ically via dust and gas in the intervening line of sight. The stel- 
lar and scattered quasar components are modeled simultaneously 
to achieve the best fit to the observed spectrum. Each model 
spectrum is masked in the same manner as the data for which it is 
intended (see §4). A model spectrum A/a is of the form: 



Ma = X, ssp.xj rx + Qx 6 (1) 

with the following definitions: 

• sspix is the stellar component of the model which consists 
of a reduced basis of 1 5 sim ple stellar populations (SSPs) from 
the iBruzual & CharlotI | |2003D , BC03, synthesis model, degraded 
to the instrument resolution, assuming solar met allicity, Padova - 
1994 models, and a Chabrier initial mass function ( IChabriej200"3h . 
Extensive simulations by CF05, have found this reduced spectral 
base (ti=0.001, 0.00316, 0.00501, 0.01, 0.02512, 0.04, 0.10152, 
0.28612, 0.64054, 0.90479, 1.434, 2.5, 5, 1 1 and 13 Gyr) to rehably 
produce observed spectral features while limiting redundancies in 
quiescent galaxies. Due to the increased complexity introduced by 
the scattered light subtraction, only solar metallicity is considered. 
Each SSP is normalized such that the area under the unmasked 
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2.0x10'* 4.0x10^ e.OxlO'* 8.0x10^ 1.0x10^ 1.2x10^ 1.4x10^ 

Number of calls 

Figure 2. Parameter space search for minimum \^ value. The three densely populated regions consist of 20 simulated annealing runs each. As the number of 
calls, or calculations of x^, advances the 'tempature' is reduced to more restrictive values. These regions are separated by 10 straight downhill searches. This 
entire procedure is performed 20 times independently to determine the model's output. 



portions of the spectrum is unity. Normalization at a single wave- 
length (Ao=4020A, as prescribed by CF05) was also investigated. 
For our sample, the alteration of normalization convention does not 
affect the results of our model more than ~la. Xi is the population 
weighting factor. 

• r\ = 10"° '*'^^ is the extinction law used to model the in 
situ reddening. A\( Av, Rv) is defined by the Galactic law of 
ICardemetaLlj 19891) assuming the mean value of Rv (3.1) in the 
diffuse ISM. A single extinction is assumed for the entire stellar 
component. 

• Qa is the on-axis observed spectrum of the quasar normalized 
such that the area under the unmasked portions of the spectrum is 
unity. = ai + a2X + asX^ + a4X'^ is the scattering efficiency 
curve, which modifies the observed quasar spectrum with the goal 
of modeling the scattered quasar light observed off axis. We have 
shown that this model component is necessary to reliably remove 
scattered quasar light ( Sheinis 2002). The model assumes no host 
galaxy contamination to the on-axis spectrum. The dominance of 
the high luminosity quasar (Mv < —23) on-axis makes this as- 
sumption acceptable. 

Finally, the individual pieces are combined as prescribed by 
equation (1) and the resulting model spectrum, A/a, is normalized, 
dividing by a single constant, such that the area under the 
unmasked portions of the spectrum is unity (the normalized model 
is denoted A/An)- 



3.1.1 Parameter range 

Our model contains 20 free parameters: (1) an extinction parame- 
ter Av, (15) SSP weighting factors Xi, (4) the scattering efficiency 



curve coefficients, a[i.4]. The range of values to explore with the 
model is constrained in the case of Xi to be 0.0 to 1.0, where 
Exi = 1.0. The range of values for Av and ay^.^ were found in an 
iterative fashion. First the model was applied to th e entire datase t 
with parameter ranges infor med from the work of ISheinii ^20021) 
and Miller & Sheini j j2003h . The range was then increased as nec- 
essary to ensure the best fit was not determined by the search- 
able parameter space. For example, the highest model Av output 
obtained is 0.4; therefore, the model samples Av from 0.0 to 1.0. 
We have tested the sensitivity of the model output to these bound- 
ary conditions by doubling the range and re-running the model on 
the observed data. No significant change in results is found. 



3.1.2 Parameter search 

Having constructed a model, we then devised a reliable means to 
find the minimum chi-squared fit. The fitting function is defined 
by: 



2 ^ / OAn - A/a„ 



(2) 



Oxn is the observed off-axis spectrum, normalized in the same 
manner as A/An, and cta is the observed spectral noise. The noise 
spectrum is generated by measuring the standard deviation of OAn 
at three to four locations devoid of prominent spectral features. 
These measurements are then extrapolated across the wavelength 
range of OAn by a low order fitting polynomial. 

The starting point in the parameter search is defined in the 
following manner. The initial stellar component is set to be 100% 
5Gyr SSP. The reddening is randomly selected, constrained by its 
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50 100 50 100 50 100 6 7 8 9 10 -0.2 0.2 0.5 1.0 

Input [%] Input x, [%] Input Xq [%] Input S x, log t; Input B-V 



Figure 3. Monte-Carlo simulation results for ideal noiseless off-axis spectra with 60% scattered quasar light contamination. 65 synthetic spectra, S\, are 
analysed by the model to determine the reliability of the results. For each S\, the known young, intermediate and old population percentage and the known 
{log ti,) and rest frame B-V input is compared to the model output as represented by 65 plus signs in each box. Ideally, all points would fall on the diagonal, 
denoted by a dashed line. The mean offset from the ideal and the standard deviation of the model output are displayed numerically in the upper left comer and 
graphically by the blue strip. 



range. The scattering coefficients are also randomly selected but 
then scaled so that scattered quasar hght comprises 60% of the 
off-axis model spectrum. Parameters a[i.4] are further constrained 
by rejecting any solutions which drive the resulting off-axis quasar 
flux, Q\ , below zero. 

Simulated annealing optimization jPressH 19921) is utilized to 
search the 20-D parameter space. Simulated annealing evaluates the 
merit function, in this case, and travels downhill in merit space 
to find a minimum within predefined tolerances. To help avoid lo- 
cal minima, a user defined parameter, traditionally called 'tempera- 
ture', is set to determine the probability of travelling uphill in merit 
space. The higher the temperature the more probability of traveling 
uphill in the merit space. A zero value for temperature corresponds 
to a straight downhill search. 

The utilization of this optimization routine, described below, 
is somewhat of an art. The temperature value must be decreased 
slowly so that the merit space can be adequately sampled, yet fast 
enough to arrive at a timely solution. After extensive experimen- 
tation the following implementation has been found to produce fa- 
vorable results. First, twenty iterations of simulated annealing are 
performed with linearly decreasing temperature. If an iteration fails 
to find a minimum value within a factor of two of the overall 
minimum then the next iteration is reset to the overall minimum; 
otherwise the next iteration starts where the previous iteration left 
off. At the completion of the simulated annealing runs, the ten best 
distinct solutions are used as input for straight downhill (tempera- 
ture = 0) runs. The current global minimum is then used as the input 
to the next twenty iterations of simulated annealing. This procedure 
is carried out for a total of three loops; within each loop the tem- 
perature parameter is set to more restrictive values. An example of 
this parameter space search is shown in Fig. 2 for a typical mini- 
mization. Overall approximately 150,000 points are sampled in the 
parameter space. This entire procedure is performed 20 times in- 
dependently. The best value of these runs, which consists of 
approximately 3 x 10*' calculations of is defined as the model 
output. 



3.1.3 Quantities of interest 

Guided by the analysis of CF05, we limit the scope of our project to 
a coarse, but well recovered, description of the stellar populations. 
Like CF05, we only attempt recovery of 'young' (t < 100 Myr), 
'intermediate' (100 Myr ^ t < 1 Gyr), and 'old' (t > 1 Gyr) stel- 
lar contributions (xy, xi, and xo, respectively). Additionally, we 
report the mean light-weighted stellar age. 



JV, =15 

{logU)j^ = ^ xjogti (3) 
1=1 

and the rest frame B-V colour as derived directly from the listed 
BC03 SSP values. These quantities have been studied extensively 
in quiescent galaxy samples and to a lesser extent in AGN host 
galaxy samples. Thus, we may assess the consistency of our 
results to previous AGN host studies. Furthermore, we may 
compare these properties of quasar host galaxies to quiescent 
galaxy samples. Analysis of other quantities of interest (see CF05) 
are deferred to future studies. 

3.2 Monte-Carlo Simulations 

The error in our method is estimated via a Monte-Carlo routine. 
Synthetic spectra, S\, are generated with known parameters and 
then degraded by noise. The same procedure used for the observed 
data is applied to determine the agreement between the known in- 
put and model output. In other words, S\n replaces 0\n in the 
fitting routine. S\ are constructed in the same manner as AI\ with 
the exception of metallicity, SSP base and added noise. Because the 
observed galaxies will likely include non-solar metallicities, even 
though the model grid does not, the metallicity is allowed to devi- 
ate from solar, and is randomly selected from Z — 0.2, 1 and 2.5 
Zq. The SSP base is extended to include all 221 BC03 instanta- 
neous bursts assuming a Chabrier initial mass function. A synthetic 
spectrum is of the form: 

/3x221 \ 

Sa = I ^ xj sspjx 1 rx + Qx Ca (4) 

Poison noise is added until the synthetic spectra obtain the cor- 
responding observed noise as measured in the 5500-5700A spec- 
tral window. The 5500-5700A spectral range was selected because 
there are no prominent emission or absorption features; addition- 
ally, all observed spectra include this range. The value of Av and 
a[i.4] are randomly selected, constrained by the applicable range. 
Parameters a[i.4] are further constrained by rejecting any solu- 
tions which drive the resulting off-axis quasar flux, Qx^\, below 
zero. The scattered quasar flux is then scaled to the contamina- 
tion percentage of interest. The value of xj is randomly selected 
to uniformly sampl^fleach stellar age bin (young, intermediate, and 

To uniformly sample each stellar age bin the following procedure is em- 
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4000 5000 6000 7000 

^0 ( A) 

Figure 4. Model fit to the 3C 273 4N pointing. Tlie observed off-axis rest frame spectrum is highlighted in cyan, with masked regions indicated in black and by 
vertical lines extending to zero flux; masked regions that are narrow, [OIII] at 5()07A for example, are most easily seen by looking at the residuals at the bottom 
of the plot. The model output is over-plotted in red. The residual is shown with diamonds. Model components are also shown. The scattered light component is 
indicated by the thin black line and comprises 68.8% of the total observed flux. The stellar component is also shown, broken up into young (blue), intermediate 
(green), and old (magenta) bins. In this example, only the old stellar component, which comprises 86.2% of the stellar flux, is clearly visible. Results for all 
other pointings are found in Appendix A, Fig. Al. 



old). Thus, synthetic off-axis spectra are generated with generalized 
metallicity and SSP base. Sixty-five unique S\ are constructed to 
test the ability of the model to recover known input parameters. We 
will use this Monte-Carlo routine first to test the associated uncer- 
tainties in a noiseless test case (§3.2. 1) and the reliability to recover 
the scattered nuclear light (§3.2.2), before we will apply it to our ac- 
tual data to determine the errors and biases in the recovered model 
parameters (§4). 



3.2. 1 Model performance given noiseless spectra 

A concern regarding the use of our multi-component model is in- 
herent degeneracies. For instance, can a young stellar population 
mimic scattered quasar light? The legitimacy of the model relies on 
the subtle differences between the various components. The red- 
dening function, which only operates on the stellar component, 
should be distinct from the scattering function which operates on 
the observed on-axis spectrum. The absorption features in the stel- 
lar component generally do not match quasar absorption lines. The 
scattering function, applied to the quasar spectrum, uniquely has 



ployed: The first of 65 synthetic spectra, Sx, is randomly assigned a xq 
value between 0-100%. xj and xy are then randomly assigned the remain- 
ing flux. The second (third) Sx is conshucted similarly, except xj (xy) 
is first assigned a value between 0-100%. This cycle is repeated beginning 
with the forth Sx- To further illustrate the procedure let us consider the 
construction of the first Sx ■ Suppose xq was assigned a value of 72%. The 
program populates xq by randomly selecting a SSP and a metallicity and 
then randomly assigning a percentage of the total xq value. This process 
continues until a:o=72% to within ±0.005%. xj andiy are then populated 
in the same manner. Thus, the stellar components of the synthetic spectra 
are constructed. 



broad-line components. However, can quantities of interest be reli- 
ably recovered given the many possible interactions between model 
components? 

To test the ideal performance of the model we con- 
ducted noiseless Monte-Carlo simulations. For this purpose, the 
Sloan Digital Sky Surv ey (SDSS) composite quasar spectrum 
jVanden Berk et al. 2001) was utilized as the on-axis quasar spec- 
trum, Qx. The scattered quasar light was scaled to a value of 60%. 
The results of this test are shown in Fig. 3. Referring to the figure, 
input versus output values of xy, xi, xq, {log t-,,) and B-V are 
shown for 65 synthetic spectra. Ideally, all points would fall on the 
diagonal, denoted by a dashed line (where input equals output). The 
mean offset from the ideal, mean(output-input), plus the standard 
deviation of the model output, ±stdev(output-input), are displayed 
numerically in the upper left corner and graphically by the blue 
strip. The fitting routine reliably recovers {log ti,)j^ and B-V in the 
ideal case of noiseless spectra. The values of xy, xj and xq are 
less well constrained, with systematic offsets of up to 6% and one 
sigma deviations of up to 13.4%. 



3.2.2 Scattered light verification tests 

Additional Monte-Carlo tests were conducted to check for consis- 
tency in the model. Inputting simulated off-axis spectra consisting 
of 100% scattered quasar light and a signal-to-noise of 10, the al- 
gorithm reliably assigned no flux to the stellar component (65 iter- 
ations resulted in a mean result of 99.92±0.01% scattered quasar 
light). Inputting 100% stellar component (S/N=10), the algorithm 
reliably assigned no flux to the scattered light component (65 it- 
erations resulted in a mean result of 2.32±1.59% scattered quasar 
light). Finally, we tested the performance of the model when con- 
strained to have a steep blue scattered quasar component, as ob- 
served in the off-axis spectrum of PHL 909. The error estimates 
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were not significantly affected by this alteration; for details see Ap- 
pendix A (PHL 909). 



3.2.3 Limitations 

Errors associated with the: 1) population synthesis model (BC03) 
(template mis-match), 2) inadequacy of the scattering efficiency 
curve, 3) non-Poisson noise and 4) the application of a single sim- 
plistic extinction law are not represented in our Monte-Carlo simu- 
lations . Despite these deficiencies, similarly constructed quiescent 
models have demonstrated consistency when compared to compet- 
ing techniques (e.g. see CF05). In §5 we address this issue by com- 
paring our results to recent studies of quasar host galaxies. Finally, 
it should be noted that the error estimates obtained reflect the one 
sigma performance of the model. Although currently applied to a 
relatively small sample, stronger error statements will be possible 
with the increase of our sample. With this in mind, when possible, 
we consider our objects as an ensemble. 



4 APPLICATION OF MODEL AND SIMULATIONS TO 
THE TEST SAMPLE 

In this section we apply our model to the test sample to constrain 
the stellar age and host colour. We then conduct 10 (one for each 
pointing) tailored Monte-Carlo simulations to estimate the reliabil- 
ity of our results. 

The following caveats should be considered in the application 
of our model to the test sample. For each object, concurrent on and 
off-axis observations have been obtained, with the exception of the 
3S pointing of 4C31.63 (At ~ 1.5yr). The possible variability of 
this quasar over this time period is not accounted for by the model. 
Line-of-sight stellar motions are not modeled for Keck objects. Our 
instrument resolution, ~ 300 km s^^, being gr eater or of the 
order of the measured velocity dispersions (Wolf & Sheinisll2008t) 
should result in minimal errors. For WIYN objects, with a„ ~ 110 
km s~^, the measured velocity dispersions are used to smooth the 
model stellar component, sspi\. Narrow emission line regions and 
sky residuals are uniquely masked for each object, giving these re- 
gions no weight in the fitting. The broad components, and in 
some cases entire broad lines, are left unmasked, though most or 
all of the H-alpha region is masked out in many objects due to 
significant sky residuals in the red end of the spectrum. Typically, 
~ 2200 spectral data points are fit for Keck objects; ~ 1600 spec- 
tral data points are fit for WIYN objects. Fig.|4]depicts the result of 
the method applied to 3C 273 4N. The results for all other off-axis 
pointings are presented in the appendix. Fig. Al. Table 3 summa- 
rizes these findings. 

The reliability of these results are then determined with tai- 
lored simulations. Synthetic spectra are generated with features 
representative of the observed spectra by (1) degrading S\ to the 
observed noise level, (2) applying the same mask constructed for 
the observed spectrum, (3) using the appropriate observed on-axis 
spectrum as Q\ and (4) scaling the scattered quasar light to the 
observed contamination level as determined by the model outpu(f[ 
An example of a simulation with the same on-axis spectrum, noise 



* The contamination percentage is well recovered as tested by the Monte- 
Carlo routines; for example, the contamination percentage for PKS 0736- 
017 is recovered within ±3.6%. This observation has the most limiting 
signal-to-noise and provides an estimate of the error upper limit. 



characteristics, applied mask and percent scattered quasar light as 
the observed off-axis spectrum of PKS 0736-1-017 is displayed in 
Fig. 5 top row. The bottom row of Fig. 5 presents the simulation 
results for PG 1309-1-355. The PKS 0736-1-017 off-axis pointing has 
the worst signal-to-noise, while the PG 1309-1-355 off-axis point- 
ing has ^average signal-to-noise and quasar contamination. The 
results of all simulations are shown in Table 3. Throughout the rest 
of the paper, mean offsets are applied to model results and Monte- 
Carlo one sigma errors are quoted. 



5 COMPARISON TO RECENT STUDIES 

We compare our spectroscopically derived rest frame B-V host 
colours to the available results from imaging studies. These imag- 
ing studies spatially subtract the point spread function of the AGN 
to reveal the host galaxy. Our method, which spectrally subtracts 
the scattered quasar light, provides a complimentary technique to 
constrain the host B-V colours. 

Our off-axis method probes a population of luminous quasars 
unattainable by most other spectroscopic techniques. The stellar 
properties of this population have yet to be determined for a large 
well defined sample. This deficiency is due to the overwhelming 
emiss ion from the central quasar. Previous off-axis studies (e.g. 
Hughe s et aH lioOO) have dealt with this problem by carefully se- 
lecting distant (~5") off-axis slit positions. This can significantly 
reduce the typical scattered light observed off axis (~20%). Our 
approach is to more rigorously model the quasar light scattered off 
axis. This allows for 1) scattered light fractions of up to ^--^80% and 
2) more freedom in slit positioning. Furthermore, our method draws 
from established quiescent population synthesis models to more ro- 
bustly model the stellar content and to estimate the reliability of our 
results. 

Recent advances in on-axis observation al techniques have 
allowed for the st udy of luminous quasars jjahnke et al.l |2007| ; 
iLetawe et al]|2007h . For objects with Mv{nuc) > My (host), 
our method can compliment these studies by observing off-axis 
(> Qkpc) stellar content with favorable signal-to-noise. We note 
that on-axis studies require nucleus-to-host ratios of < 15 to detect 
stellar absorption features (L07). Objects such as NAB 0205-1-02, 
HE 0530-3755 and 3C 273 cannot be observed on-axis due to this 
criterion. Off-axis studies do not have this limitation. Thus, our 
method is needed to obtain optical stellar spectra for these objects. 

In the following subsections we compare our study to previ- 
ous imaging and spectroscopic work in more detail. Furthermore, 
we quantify how our sample compares to applicable SDSS AGN 
studies. 



5.1 Recent photometric studies 

Host B-V c olours have b een re ported for nearby (z < 0.197) 
quasars by Ijahnke et al.l i2004h. hereafte r J04, and for nearby 
(z < 0.139) BL Lacs bv lHvvonen et alj j2007l) . hereafter H07. 
J04, studying multicolour data from a sample of AGN distributed 
around the classical dividing line in luminosity between Seyfert 1 
galaxies and quasars, found bulge dominated hosts to have bluer 
colours than expected compared to their quiescent counterparts; 
disc hosts were also on average bluer, but not by a significant 
amount. H07, also studying multicolour data, found their exclu- 
sively elliptical BL Lac hosts to be on average bluer than expected 
compared to their quiescent equivalent. Table 4 shows the mean 
B-V colours for the host galaxies in J04 and H07, as well as in 
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Pointing 



Young stell 
pop - % flux 



Inter stell 
pop - % flux 



Old stell 
pop - % flux 



B-V 



3C273 4N 
4C 31.63 2N 
4C 31.63 2.5E 
4C 31.63 3S^ 
PKS 1302-102 2.3N 
PKS 0736+017 4.5NW 
PKS 2135-147 3W 
PKS 2349-014 4N 
PKS 2349-014 3S 
PG 1309+355 4.5SW 
PHL909 4.5N 
PG 0052+251 3S 
PG 1444+407 3S 



13.8(-3.0±16.9) 
15.7 (-1.9±16.2) 
18.3(-1.1±16.4) 
0.0(-2.1±17.0) 
0.0 (-5.6±27.0) 
0.0 (+4.2±21.3) 
0.0 (-2.0±22.6) 
6.3 (-2.5±11.4) 
11.1 (-4.7±12.7) 
11.7 (-2.1±14.4) 
0.0(-8.3±18.4) 
43.3 (-8.0±19.2) 
9.5 (-7.8±23.8) 



0.0 (+0.3±19.1) 
0.0 (-1.8±16.5) 
0.5 (-2.3±17.7) 

25.3 (-1.0±18.2) 

12.4 (+4.3±26.7) 
0.0 (-5.7±23.2) 
0.0 (-0.6±23.7) 
0.8 (+0.0±11.3) 

15.3 (+0.8±13.2) 
9.5 (+1.9±16.1) 

43.8 (+4.9±20.0) 

16.5 (+4.5±20.9) 
58.2 (+2.6±27.6) 



86.2(+2.6±13.8) 
84.3 (+3.7±9.6) 
81.3 (+3.4±11.8) 
74.7 (+3.1 ±9.9) 
87.6(+1.3±14.8) 
100.0 (+1.5±12.4) 
100.0 (+2.6± 12.0) 
92.9 (+2.5±9.6) 
73.6 (+3.9±8.8) 
78.9(+0.2±10.7) 
56.2 (+3.4±12.4) 

40.1 (+3.6 ±11.5) 

32.2 (+5.2±17.7) 



9.40 (-0, 
9.00 (+0. 
9.18 (+0. 
9.73 (+0, 
9.34 (+0, 
9.56 (-0, 
9.52 (+0, 
9.52 (+0. 
9.54 (+0. 
9.33 (+0. 
9.45 (+0, 
8.09 (+0, 
8.86 (+0, 



.01 ±0.24) 
,06±0.28) 
).00±0.32) 
.03±0.26) 
.09±0.49) 
.05±0.32) 
.01±0.34) 
l.01±0.20) 
.06±0.21) 
.02±0.26) 
.14±0.45) 
.09±0.34) 
.15±0.42) 



0.77 (-0, 
0.69 (+0, 
0.74 (-0, 
0.83 (+0 
0.76 (+0 
0.83 (-0 
0.82 (-0, 
0.79 (+0, 
0.78 (+0 
0.76 (+0, 
0.73 (+0 
0.32 (+0 
0.52 (+0 



00±0.05) 
01 ±0.06) 
00±0.05) 
.01±0.05) 
,02±0.12) 
,01 ±0.07) 
00±0.07) 
00±0.04) 
,00±0.03) 
00±0.06) 
,02±0.10) 
.03±0.05) 
.03±0.09) 



Table 3. Modeled off-axis stellar populations, {log ti,) and rest frame colour are shown. Host stellar population offset and ± Icr error are shown in parentheses 
as estimated with Monte-Carlo simulations. A denotes the observation with non-concurrent on and off-axis observations. 




Figure 5. Top row: Monte-Carlo simulation for PKS 0736+017 (off-axis pointing of most limiting signal-to-noise). Bottom row: Monte-Carlo simulation 
for PG 1309+355 (off-axis pointing of ~average signal-to-noise and quasar contamination). 65 synthetic spectra are analysed by the model to determine the 
reliability of the results. For each Sx, the known young, intermediate and old population percentage and the known (log t*}^ and rest frame B-V input is 
compared to the model output as represented by 65 plus signs in each box. Ideally, all points would fall on the diagonal, denoted by a dashed line. The mean 
offset from the ideal and the standard deviation of the model output are displayed numeiically in the upper left comer and graphically by the blue strip. 



our s ample and for quiescent galaxies of va rious morphologies 
(from [pukugita. Shimasaku & Ichikawalll995h . In agreement with 
J04 and H07, we find the rest frame B-V colour for the elliptical 
host galaxies of our sample to be bluer than one would expect given 
their quiescent morphological counterparts. This provides confir- 
mation that elliptical quasar hosts are distinguishable in B-V colour 
from inactive elliptical galaxies using a spectroscopic technique. 
Spirals hosts for our sample are found to be consistent with S c 
galaxies which have a B-V colour of ~0.54 jFukugita et^ll995h . 
A comparison to quiescent galaxies is also presented vi a a colour- 
magnitude diagram, Fig. 6. Quiescent COMBO- 17 jWolf et al.l 
|2004) late-type galaxies (black '+'), and quiescent COMBO-17 
early-type galaxies (black '*') are shown in relation to J04 hosts 
(blue lower-case letters; 's' for spiral or disc, 'e' for elliptical mor- 
phology), H07 hosts (red 'e'), and this work (green upper-case let- 
ters; 'S' for spiral or disc, 'E' for elliptical morphology). Quiescent 
early and late-type galaxy designations are d etermined by the U-V 
vs Mv cut as proposed bv lBell et alj j2004h for z=0.25. A typical 
error bar is shown for our sample in the lower right of Fig. 6. 

Although we confirm bluer elliptical host galaxy colours, we 
find that our sample is on average much redder than J04. This off- 
set may be due to the different effective radii probed since our 



method is only sensitive to the stellar content observed off axis. 
For example, if nuclear starbursts were a common occurrence in 
quasar hosts, then one might expect the outer st ellar populations 
to exh ibit redder colours. However, results from lKauffmann et al.l 
1 I2OO3I) suggest that the different radii probed may not be responsi- 
ble for this observed offset. Studying a large sample of obscured 
AGN hosts, they find that star formation in the hosts of power- 
ful AGN (approximately one magnitude fainter than a typical ob- 
ject in this study) is not concentrated primarily in the nuclear re- 
gions; it is spread out over scales of at least several kilo-parsecs. 
Furthermore, the host galaxies in the H07 study exhibit a wide 
range of host colour gradients, with the majority displaying bluer 
colours off axis. T his negative colour gradient is also observed 
for radio galaxies iGovoni et al.l [2OO0I) and quiescent ellipticals 
jPeletier. Valentiin & Jamesonlll99 (i). Assuming these conclusions 
hold true for quasar hosts, one would expect our sample to be bi- 
ased toward bluer colours rather than the observed red offset. 



An alternative explanation for our redder host colours com- 
pared to the J04 study is that our sample is on average more radio 
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Figure 6. Rest frame colour-magnitude diagram for local quiescent galaxies and active host g alaxies. Black plus signs indicate (z<0.4) quiescent late-type 
galaxies, while black asterisks indicate (z<0.4) quiescent early-type galaxies iWolf et al]|2004l) . J04 quasar hosts are blue lower case letters ('s' for spiral or 
disc, 'e' for elliptical morphology). H07 BL Lac hosts are all elliptical and are indicated by red 'e\ Quasar hosts from this work are indicated by green upper 
case letters ('S' for spiral or disc, 'E' for elliptical morphology). A typical error bar is shown for our sample in the lower right. '1' and '2' denote objects with 
multiple pointings, 4C 31.63 and PKS 2349-014 respectively. A denotes the observation with non-concurrent on and off-axis observations. 



The only J04 object which overlaps in B-V colour with our 
main sample is their most radio loud. The two objects from our 
sample which overlap in colour with the main J04 sample are also 
our most radio quiet. This relation of colour to radio luminosity is 
shown in Fig. 7. Above L^ohz ~ lO'*" erg s~^, the host galax- 
ies are systematically redder. Given the correlation between radio 
luminosity and bulge mass found for our sample (Wolf & Sheinis 



2008, LsGHz 



M, 



3.56 
bulge 



), L^GHz = 10''^ erg s ^ corresponds 



to Mbuige ~ 10 Mq. The rms scatter of the correlation is 1.09 
dex. This trend f or radio lou d objects to be located in redder hosts 
is not surprising: iBest et"al] jlOOS) found that radio-loud AGN are 
preferentially located in older more massive galaxies. However the 
distinct transition we see in our sample should be investigated on a 
larger sample. 



Sample 



B-V 



E Quasars (this work) 
S Quasars (this work) 
E Quasars (J04) 
S Quasars (J04) 



0.78±0.01 
0.55±0.12 
0.52±0.05 
0.56±0.01 



Quasars LsoHz > 10*° erg s"i(this work plus J04) 0.77±0.02 

Quasars LsgHz < 10"*" erg s"^(this work plus J04) 0.50±0.02 
E BL Lacs (H07) 0.8±0.2 
E (Fukugita et al. 1995) 0.96 
Sb-Sc (Fuku gita et al. 1995) 0.57 
Sc-Sd (Fukug ita et al. 1995) 0.50 

Table 4. Available rest frame B-V host mean colours and quiescent mean 
colours. Reported eiTors are standard deviations of the mean. 'E' denotes 
elliptical and 'S' denotes spiral. 



How can this be understood in terms of stellar ages? The 
m ean light- weighted age is correlated to various colours (e.g., u- 
r, iMateus et alfcooeh . We verify this correlation for our colour of 

^ Ra dio luminosity is estimated utilizing NED and NVSS iCondon et alj 
Il998l) . A spectral index of -0.5 is assumed. Upper limits are based on the 
sensitivity of NVSS ~ 2.3 mjy @ L4 GHz. 



interest, B-V, by constructing 65 test spectra, Fig. 8. Using this re- 
lation we convert the B-V colours of J04 to mean light- weighted 
ages. Fig. 9. While the observed colour trend is more pronounced, 
it is found that all {log i*)^ values above L^ghz =10'*" erg 
are older than valu es below this threshold. Previous studies (CF05; 
iMateus et d][200"6) have studied the mean light-weighted age and 
found this quantity to be frequently associated with the most 
recent epoch of starbursts. Assuming major merger progenitors 
for our quasar hos ts with a significant past starburst (e.g., see 
iHopkins et al.ll20()6l) , {log ti,) may indicate the approximate tim- 
ing of the merger event. It should be noted that our definition of 
{log tt,) J dif f ers sli ghtly from the definition adopted by CF05 and 
IMateus etal] ( |2006|) due to different SSP normalization conven- 
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Figure 7. Radio luminosity as a function of rest frame B-V colour. Green 
upper case letters indicate the results of this study, blue lower case letters 
indicate the results of J04. '1' and '2' shown above denote objects with 
multiple pointings, 4C 31.63 and PKS 2349-014 respectively. A denotes 
the observation with non-concurrent on and off-axis observations. H07's 
BL Lacs are not shown due to the expected significant beaming of radio 
emission. 



Figure 9. Radio luminosity as a function of mean stellar age weighted by 
flux. Green upper case letters indicate the results of this study, blue lower 
case letters indicate the results of J04. '1' and '2' shown above denote ob- 
jects with multiple pointings, 4C 31.63 and PKS 2349-014 respectively. A 
denotes the observation with non-concurrent on and off-axis observations. 



radio bright sample and (t*)^ = 
sample. 



570 
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Figure 8. Relation between mean light-weighted stellar age and rest frame 
B-V colour derived from 65 randomly constructed stellar components, 

Els 
Xi sspix- 



tions (see §3.1). For our sample, conversion to this alternative defi- 
nition modifies our reported {log ti,) results by ~1(t or less. 

To estimate this time-scale we divide the combined sample 
into two groups; a radio bright sample with Lsghz > 10'"' erg s~ ^ 
displaying an average B-V colour of 0.77 and a radio faint sample 
with L^GHz < lO'^'^erg displaying an average B-V colour of 
0.50 (average colour for this grouping is also reported in Table 4). 
LzjGHz is found to correspond to bulge mass and may provide a 
more significant division than morphology whose classification is 
hampered by quasar point source subtraction and un-modeled tidal 
features. We find that, on average, (t*)^ — 2.7^0 5 Gyr^for the 



5.2 Recent spectroscopic study with imaging 

This result for r adio loud object s comp ares remarkably we l l to re - 
cent findings o f Canalizo et al. I bOOa) and iBennert et al] ( l2008h . 
ICanaUzo et^ ( l2006h report preliminary deep Keck LRIS spectral 
analysis of luminous z ~ 0.2 quasar hosts which finds evidence for 
major starburst episodes with ages ranging from 0.6 Gyr to 2.2 Gyr. 
Incidentally at least two of the host galaxies studied by Canalizo & 
Stockton (in preparation) overlap with our sample, PKS 0736-1-017 
and PHL 909. They find that both of these objects located in our 
radio bright group have massive 2.2 Gyr starbursts, consistent with 
our colour and mean light- weighted age estimate. However, see Ap- 
pendix A for an inconsistency in the derived stellar populations for 
PHL 909. lBennert et alj j2008l) . again studying PKS 0736-1-017 and 
PHL 909 along with three other hosts in deep HST ACS images, 
find shells and tidal tails indicative of merger events. Consulting 
N-body simulations they conclude that in general the observed fine 
structure can be explained by either a recent minor merger or an 
older major merger '--^IGyr ago. 

Implicati ons of a ^Gyr posts tarbu rst quasar are discu ssed 
extensively in ICanalizo et al j j2007h and iBeimert et al. I ( l2008h . In 
summary the curren tly estimated (opticall y luminous) quasar life- 
time (10*^ - lO^yrs jYu & Tremainei2002h ) appears to be too short 
to be triggered at the same time as the starburst. iBennert et al.l 
(i2008J) conclude that this may imply a scenario in which quasar 
activity is significantly delayed or a scenario in which intermit- 
tent quasar activity encompasses a longer duration. Alternatively, a 
more recent minor merger could explain the observed quasar. Our 
results do not rule out this scenario. However, the average mass-to- 
light ratio for our sample being ~10 l lWolf & SheinisH2008) . con- 
strains any recent (t^lGyr) starburst to be relatively insignificant 
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5.3 Off-axis studies 



Sample 



{log 



The off-axis (5") spectroscopic study of'Hugh es et alj ( l2000l here- 
after HOO) consists of three matched subsamples of nine radio quiet 
quasars (RQQs), ten radio loucQ quasars (RLQs), and seven ra- 
dio galaxies (RGs). Similar to our study, all quasars are luminous 
{Mv < -23) and local (0.1 < z < 0.3). 

HOO determine the 4000 A break strength, an indicator of 
the mean stellar age, for each object. The average values for RQQ 
hosts, RLQ hosts, and RGs all indicate younger mean stellar ages 
than local inactive elliptical galaxies, in agreement with our study. 
However, HOO notes wide scatter in their measured values. Addi- 
tionally, measurements indicative of younger ages were typically 
associated with either poor signal-to-noise data or with spectra 
that show significant scattered quasar light. Thus HOO defer fur- 
ther analysis of the stellar populati ons to the followup spectro- 
photometric modeling conducted by iNolan et al. (2001, hereafter 
NOl). 

NOl develop two population synthesis models. The first model 
consists of a 0. 1 Gyr SSP and a second SSP component whose age 
and contribution are allowed to vary to obtain the best fit. Like 
our study, solar metallicity is assumed. However, two SSPs are used 
to construct the stellar population in contrast to our model which 
incorporates a base of 15 BC03 SSPs. NOl also develop a three 
component model which uses the observed nuclear spectrum of the 
RQQ, 0054-1-144, to account for scattered quasar light detected off 
axis. While our approach also uses observed nuclear light to model 
off-axis scattered light, our model does not assume a universal nu- 
clear spectrum. The observed nuclear spectrum is matched to the 
off-axis observation (e.g. the model input for the quasar contamina- 
tion in the 3C 273 off-axis spectrum is the observed 3C 273 nuclear 
spectrum). Additionally, it should be noted that whenever possible 
(9 out of 10 pointings) we have observed on and off-axis obser- 
vations simultaneously to avoid variation in the quasar spectrum 
which could hamper an accurate scattered-light subtraction. As a 
final point, we allow for the continuum shape of the nuclear light 
observed off axis to be altered by the scattering efficiency curve, 
which we have show n is necessary for reliable scattered light re- 
moval ( ISheinisll2002l) . 

Interpretation of the NOl model output indicates that the stel- 
lar mass of all three subsamples studied (RQQ, RLQ, and RG) is 
dominated by old stars of age 8-14 Gyr. This result is used to sup- 
port the claim that quasar hosts are "to the first order, indistinguish- 
able from 'normal' quiescent giant elliptical galaxies." If we take 
NOl's derived age to be representative of the mean mass-weighted 
stellar age, then it is not clear that this quantit y can be used to differ - 
entiate between late and early type galaxies. iMateus et ai] l l2006h . 
studying 50,000 luminous galaxies from the SDSS, show that the 
distribution of mean mass-weighted stellar ages is not bimodal, as 
both star-forming and passive galaxies have formed a large frac- 
tion of their stellar mass at early times (~ 12.4 Gyr in the past). 
We may comment on the result of NOl by referring to the mass- 
to-light ratio derived fr om the measured stellar velocity dispersion. 
IWolf & Sheinisl ( l2008h found that the average mass-to-light ratio 
for our quasar host sample is ~ 10, which supports the domi- 
nance by mass of old stars. However, as stated in §5.1 our colour 
/ (logti,)^ determination indicates that quasar hosts are bluer / 
younger than quiescent ellipticals. 



Quasars L^qhz > W^^erg (this work plus J04) 
Quasars L^qu^ < 10'"' erg (this work plus J04) 
Passive (M06) 
AGN hosts (M06) 
'green valley' (M06) 



Star-forming (M06) 



9.4±0.1 
8.8±0.1 
9.86 
9.66 
9.53 
8.91 



Radio loud in HOO is defined as LsGjY J. > lO'^^W Hz'^ sr- 
is roughly consistent with our adopted definition. 



, which 



Table 5. Comparison of derived {log t^) values to those found in the 
SDSS study conducted by M06. Reported errors are standard deviations of 
the mean. 



5.4 On-axis studies 

J07 present a method that spatially subtracts the quasar point spread 
function from the on-axis spectrum to extract the host galaxy light. 
Unlike our method, which is constrained to collect light > 9kpc off 
axis, this technique allows for the inner regions of the host galaxy 
to be probed. Like our study, the J07 test sample consists of nearby 
(z < 0.3) quasars. However, this sample consists of quasars that are 
typically less luminous in optical and radio bands than our study. 
J07 model 13 host spectra with a two component BC03 SSP model. 
Eight of the spectra were deemed trustworthy, yielding a typical 
light-weighted stellar age of 1-2 Gyr. For a rough comparison, our 
sample average (i*)^ is approximately 2 Gyr. If we limit our sam- 
ple to radio quiet objects, to more closely match the sample of J07, 
we find an average (f*)^ of 1 Gyr. Confidence in this comparison 
is hampered by the small number of objects probed and the sample 
mismatch: our targets are typically more luminous. 

Despite this mismatch, there is one object that is modeled by 
both groups, the radio loud quasar PKS 1302-102. While 107 find 
a stellar age < lOOMyr, our results indicate a predominantly old, 
> IGyr stellar population. This disagreement may indicate a radial 
dependence on the stellar populations or model systematics. 

L07, utilizing an on-axis method complementary to the J07 
technique, examine the host galaxies of 20 luminous quasars whose 
redshifts and nuclear optical luminosities are well matched to our 
sample. However, their optically selected sample contains far fewer 
RLQs. L07 measure diagnostic absorption and emission lines and 
then compare these to known quiescent galaxy values. Half of 
their sample was found to have young Sc-like stellar populations. 
The sample examined by our study indicates a redder and hence 
older stellar population than a typical Sc galaxy. Our studies are 
again brought into closer agreement if we consider only our RQQ 
hosts. However, no secure agreement can be claimed due to our 
small RQQ sample. Further complicating the comparison, L07 also 
model PKS 1302-102 and find evidence for a young Sc-like stellar 
population. Additional interpretation of this disagreement is dis- 
cussed in Appendix A. 



5.5 Sloan Digital Sliy Survey studies 

The Sloan Digital Sky Survey (SDSS) has been utilized to 
study AGN hosts by a number of groups. These studies ben- 
efit from a large we ll defined sample. N arrow-line AGN stud- 
ies teauffm ann et al.ir2003,. hereafter K03; IZakamska et"d]|2003l : 
iMateus etal . 2006, he reafter M06) and the broad-line AGN study 
bv lVanden BerketalJ j2006l hereafter V06) are compared to this 
work below. 

K03 study the host galaxies of 22,623 local type 2 AGN. Mod- 
ulo deficiencies in AGN unification models, quasar hosts can be 
studied via their narrow-line counter-parts. From 4000 A break 
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measurements, K03 find that high luminosity AGN have younger 
mean stellar ages than normal early-type galaxies. Additionally, 
they find a large fraction of powerful type 2 AGN have experi- 
enced significant starbursts in the past 1-2 Gyr. This is similar 
to our findings: {t*) ^ — 2.7j;Qg Gyrs for the radio bright sam- 
ple and = 570^110 Myrs for radio faint sample. However, 
there is little overlap in luminosity between our samples. The K03 
powerful (L[OIII]- 10*^ erg s"^ ) AGN are about one magnitude 
fainter than our sample's nuclear luminosities. A smaller more lu- 
minous sample of SP SS narrow-line AGN has been examined by 
IZakamskaetaLl l l2003h . In qualitative agreement with K03, they too 
find evidence for relatively blue / young host galaxies. 

M06 divide a sample of 50,000 SDSS galaxies into star- 
forming, passive, and AGN hosts spectral classes. The AGN probed 
by this study are comparable in luminosity to the K03 sample. M06 
find a clear bimodal distribution for the mean light-weighted stel- 
lar age among star-forming and passive galaxies. As summarized in 
Table 5, median values of {log t-i) ^ = 8.91 for star-forming galax- 
ies and {log t-,,) = 9.86 for passive galaxies are found. The green 
valley is located at {logtt,)^ ~ 9.53 and the median {logti,)j^ 
for AGN hosts is 9.66. The M06 median {log ti,) for AGN hosts 
is heavily weighted toward the more numerous low luminosity ob- 
jects. In agreement with K03, they find that the typical {log t*)^ 
decreases for more luminous AGN. Our radio bright sample with 
{logt-),)^ ~ 9.4 ± 0.1 is in close proximity to the green valley 
(slightly younger and bluer). We find the radio faint sample with 
{log ti,) ~ 8.8 ±0.1 to fall very close to the peak of star-forming 
distribution. Qualitatively we find this to be in overall agreement 
with our colour-magnitude diagram. Fig. 6. 

V06 study the host galaxies of 4,666 local broad-line AGN. 
Like K03 and M06, they find that higher luminosity AGN hosts are 
bluer /younger than normal early-type galaxies. This sample con- 
tains ~ 1750 quasars with luminosities comparable to our study. 
However, the model employed by V06 fails when the on-axis host 
flux fraction falls below 10% or when the signal-to-noise is less 
than 10. These criteria prevent the study of ~ 70% of their lumi- 
nous objects and would exclude the majority of the objects (7 of 
10) in this study. 



6 SUMMARY AND CONCLUSIONS 

We have presented an off-axis technique to spectroscopically con- 
strain the colour and the stellar ages of quasar {Mv{nuc) < —23) 
host galaxies. Our method draws heavily from the quiescent galaxy 
model of CF05, utilizing a basis of BC03 SSPs and Cardelli et ah 
( Il989l) dust extinction in a similar fashion. Complicating our fit- 
ting routine is the residual scattered quasar light (a combination of 
atmospheric, instrumental and host galaxy scattered light) which 
must be accounted for in the off-axis spectrum. Scattered light is 
modeled by altering the the observed nuclear spectrum by a low or- 
der polynomial, while simultaneously fitting the constituent stellar 
populations of the host galaxy. Furthermore, Monte-Carlo simula- 
tions are tailored to each observed pointing. Thus, the ability of the 
model to recover known parameters from synthetic spectra is de- 
termined. It is found that quasar host B-V colour and {log tt)j^ are 
well recovered. These parameters are then compared to previous 
studies. Overall consensus is found giving further credence to our 
model (however, see Appendix A for discussion of individual cases 
that potentially disagree with previous results). 

Our method probes a population of luminous quasars unattain- 
able by most other techniques. The stellar properties of this popu- 



lation have yet to be determined for a large well defined sample. 
For the smaller sample (Nobj ~ 20) studies that can access this 
luminous regime, our method can be used in conjunction to gain 
further insight. For example, our technique can observe the off-axis 
(Roba^ 9kpc) stellar content of objects (nucleus-to-host ratio>15) 
inaccessible to the on-axis studies of 107 and L07. 

The full potential of this method has yet to be realized. To 
date, we have analyzed the properties of 10 quasar host galaxies. 
With this sample, we confirm that elliptical quasar hosts are distin- 
guishable (bluer) from inactive ellipticals in rest frame B-V colour. 
Additionally, we note a trend for radio luminous quasars to be lo- 
cated in redder host galaxies in comparison to their less radio lu- 
minous counterparts. However, general conclusions await further 
observations. Currently we have 10 additional WIYN off-axis ob- 
servations in the process of being reduced. Furthermore, time has 
been allotted on the WIYN telescope to increase our total analyzed 
sample size to ~ 25 objects. This sample will consist of roughly 
an equal number of RLQs and RQQs and allow for 1) more rigor- 
ous comparison to other studies due to the increase in sample size, 

2) further verification of our technique with the targeting of L07 
overlap objects (specifically HE 0914-0031 and HE 0956-0720), 

3) investigation of the trends found in our preliminary sample. In 
the near future, the Robert Stobie Spectrograph (RSS) on the 1 1-m 
Southern African Large Telescope (SALT) will be employed in our 
ongoing observational campaign. This phase of the project will in- 
crease the observable sample by its location in the Southern hemi- 
sphere and by allowing for higher redshift objects to be targeted 
(z ~ 0.5). 

Future work will build on the current model to analyze the 
narrow-line emission. This will provide insight into the star forma- 
tion rate, the emission mechanism, the metallicity, and the stellar 
content (e.g. see L07). The accomplishment of this program will 
provide important observational constraints on luminous quasars. 
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APPENDIX A: COMMENTS ON INDIVIDUAL HOSTS 

In this section, we describe the stellar populations found for in- 
dividual hosts and compare these to previous studies. Objects are 
listed in order of greatest 5 GHz luminosity to least, as in Table 
1. Unless otherwise mentioned: young = (t < 100 Myr), interme- 
diate = (100 Myr ^ t ^ 1 Gyr), and old = (t > 1 Gyr), where t 
is the age of the starburst. Additionally, stellar populations (young, 
intermediate, and old) are reported by percent flux contribution. 
Results from the Monte-Carlo simulations, namely the mean offset 
and the one sigma errors, are utilized in the following comments. 
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By convention, application of offsets may not force the percent flux 
contribution to be negative or greater than 100%. For output param- 
eters without offsets applied see Table 3. Fig. 4 and Fig. Al are also 
shown without offsets applied. As mentioned in §3.2.3, error esti- 
mates reflect one sigma performance; thus we expect statistically 4 
of our 13 pointings to lie beyond our error limits. 

3C 273 (PG 1226+023) - The host gala xy is classi- 
fied a s an elliptical galax y via HST imaging jSahcall et al.l 
1997; Hamilton etalj |2002|) . The off- axis spectral analysis of 
Boroson et al.. (19851) qualitatively found evidence for a significant 
contribution from an old stellar population. Our results, 88. 8± 13.8 
% flux contribution from an old stellar population at an observed 
radius of 1 1.8 kpc, are consistent with these studies. Somewhat be- 
low the 1(7 error, a young stellar component is found by the routine. 
Residual Oil and 0111 narrow line emission is noted. 

4C 31.63 (Q2201+315) - 4C 31.63 is a radio loud quasar 
whose host galaxy has been classified as an elliptical by HST imag- 
ing (Bahcall et al. 1997; Hamilton et al. 2002) and by near-infrared 
adaptive optics imaging l lGuvon et al.ll2006l) . Our results for the 
Northern and Eastern pointings are consistent with a predominantly 
old stellar population. Toward the south we find evidence for an in- 
termediate stellar population not found in the other pointings. Note 
that the 3S pointing has non-concurrent on and off-axis observa- 
tions, see §4 for details. No narrow line emission is found for this 
object. 

PKS 1302-102 (HE 1302-1017) - The host galaxy of this 
radio loud quasar has been classified as a disturbed elliptical 
by various studies jHutchings & Nefflll992l; iBahcall et alj|l997l ; 
Hamilton et"ar]|2002| - lGuvon et alj|2006l) . Ijahnke et al] ( l2007b and 



Letawe et al 



1 20071) have found evidence from o n-axis spec- 
trosco py for a young spiral-like stellar population. Ijahnke et al.l 
( 12007 ^ note that there are almost no stellar absorption lines with 
the exception of a weak Ca // K and Mg / line. This is not the case 
in our observed off-axis spectrum; Ca // K A3933 and H A3969 
are prominent with clear Mg I b A5175 and Na D A5895 absorp- 
tion features. The stellar component is found to be best fit by a 
88.9±14.8% old po pulation with a stati stical ly insignificant inter - 
mediate population. Ijahnke et al.l ( l2007h and iLetawe et al.l j2007h 
analyse the stellar content with on-axis spectra, potentially prob- 
ing much closer in to the central SMBH. Thus a plausible expla- 
nation for this disagreement is the different effective radii probed 
by our studies; however, other systematics are not ruled out. In our 
observed off-axis spectrum, scattered light from the central quasar 
is significant ~68.4%, effectively decreasing the signal-to-noise 
of the stellar component. Additionally, the B aimer decrement has 
been masked due to nebular emission, potentially removing im- 
portant spectral information. On the other hand, we estimate the 
nucleus-to-host ratio to be 13.5, which is relatively close to the 
stated upper limit of ~ 15 to reliably detect absorption features 
with the on-axis method (L07). Regardless, the observed prominent 
Ca H and K absorption lines support the claim that a significant old 
population is present at the radius probed ~ 13 Kpc. Residual Oil 
and OIII narrow line emission is noted. 

PKS 0736+017 - The host galaxy has been classified as el- 
liptical tW right, McHardv & Abraham 1998 McL ure et al.lll999l ; 
iFalomo & Ulrich ,2000, ; Hamilton et al. .2002). Others describe 
the host as a distur bed elliptical based on near-infrared imagin g 
( lDunlopetaI]|l993h and deep HST imaging teenn ert et ij]|2008h . 
Off-axis spectroscopic studies have found a stellar co mponent in- 
dicati ve of a 12 Gyr population (Hughe s et al.ll2000l ; iNolan et al.l 
l200lh . Canalizo & Stockton (in prep aration) find a 2.2 Gyr star- 
burst in this host teennert et alj|200a) which falls in our 'old' age 



bin. Our results indicate 100± 12.4% flux contribution due to an old 
population. No narrow line emission is found for this object. 

PKS 2135-147 - The host galax y of this radio loud q uasar 
has been classified as elliptical by iBahcall et alj ( I1997 ') and 
iHamiltone t al. ( 2002). Our results are consistent with the expected 
stellar population of an elliptical galaxy. We find 100±12.0% old 
stellar population at a radius of approximately 12.5 Kpc. Multiple 
residual narrow emission lines are noted for this object. 

PKS 2349-014 - The host galaxy of this radio loud quasar has 
been cl assified a s a highly distu rbed elliptica l by many imaging 
studies teahcall et al. 1997; Gu von et alj|2006l) . An off-axis spec- 
troscopic study bv .Nolan et al., (120011) found a dominant old stellar 
populat ion of 12 Gyr with a low level young popula tion. Off-axis 
studies ( Miller & Sheinisl2003l ; rWolf & Sheinisl2008l) , utilizing the 
same data analysed in this paper, have also found evidence for a 
predominately old stellar population with a lesser young popula- 
tion. Our results are consistent with a predominantly old population 
by flux. The northern pointing (14. 12 kpc), probing large tidal arms, 
is consistent with an entirely old stellar population; whereas the 
closer-in Southern pointing (9.47 kpc) we find evidence for a small 
intermediate population (16.1±13.2%). Multiple residual narrow 
emission lines are noted for this object. 

PG 1309+355 - HST imaging studies classify this object as 
a spiral ('Bahcall et al. 1997; Hamilton et al. 2002), while the near- 
infrared study of Guvon et al, (2006.) classifies the host as an elon- 
gated elliptical. If this object is truly a spiral, then this is a very rare 
object displaying a large radio luminosity, just below our radio-loud 
threshold. Our results indicate a predominantly old population. No 
significant narrow emission lines are noted. 

PHL 909 (0054+144) - The ho st is an ellipt i cal g alaxy 
as classified by the i maging studies of IBahcall et all jl997h and 
iHamilton e t alj ( I2OO2I ) . The near-infrared imaging of Dvmlop et al.l 



(1993) detects exten ded emission towar d a companion. Recent 
deep HST imaging teennert et ^ I2OO8I) detects tidal tails and 
sheUs. iBarthell ( |2006|) . analysing the far-infrared / radio correla- 
tion, finds evidence of recent or on going star formation. Off- 
axis spectr oscopic studies hav e either classified the stellar popula- 
tion as old 1 Nolan et al]|200ll) or have reported ambiguous results 
dHughes etal.ll2000h . Canalizo & Stockton (in prep.), analysing 
~2 hr exposure Keck LRIS data, find a 12 Gyr p opulation com- 
bined with a 2.2 Gyr starburst (as reported in ICanalizo et al.l 
I2OO6I) . This would correspond to a 100% old population given our 
adopted age bins; however our results indicate 59.6±12.4% old and 
48.7±20.0% intermediate by flux. This pointing has the highest 
scattered quasar component at 77.8% by flux. The scattered light 
contamination is highest at the blue end of the spectrum where the 
4OOOA break and various stellar absorption features are located. 
Thus our Monte-Carlo error simulations, which normally are not 
constrained to have steep blue scattering efficiency curves, may 
be under estimating the actual error. To test this scenario, we per- 
formed Monte-Carlo simulations with input parameters indicative 
of PHL 909 and constrained to have a steep scattering efficiency 
curve, as found by the best fit model output. The error estimates 
and offsets were not significantly affected by this alteratior[fl If 
our observations are co-spatial with those of Canalizo & Stockton, 

* As listed in Table 1 and Table 2, the original error simulations found the 
following offsets and ilcr errors: -8.3± 18.4% for young, +4.9±20.0% 
for intermediate, +3.4 ± 12.4% for old, +0.02 ± 0. 10 for B-V colour, and 
+0.14 ± 0.45 for {log t*) i^- The steep blue scattering simulations found: 
-3.1 ± 22.0% for young, -1.8 ± 22.5% for intermediate, +4.9 ± 11.9% 
for old, +0.01 ± 0.08 for B-V colour, and +0.03 ± 0.45 for {log t^,) ^. 
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Figure Al. Model fit to the observed off-axis spectra. The meanings of the symbols and colours are the same as in Fig. 4. 
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then the difference in the results may simply be a 3ct outlier arising 
from separate observations and different analysis methods. Regard- 
less, construct ion of a solar BC03 st ellar spectrum consistent with 
the findings o f ICanaUzo et"aD bOOeh results in values for the mean 
light- weighted age andS-V colour that are consistent (within la) of 
our findings. Therefore, we believe that the main discussion regard- 
ing the more reliably recovered parameters, colour and {log 
are robust. 

PG 0052+251 - The host galaxy associated with this radio 
quiet quasar is a spiral as classified by various imaging studies 
teahcall et ahll 19971 : [Hamilton et al.ll2002l) . Our results are consis- 
tent with a significant young population of 35. 3± 19.2% by flux and 
an old population of 43.7±11.5%. Multiple residual narrow emis- 
sion lines are noted for this object. 

PG 1444+4 07 - The host of this r a dio quiet quasar is c lassi- 
fied as a spiral bv lHamilton et alj j2002l) . lBahcall et all fl997h note 
that the host has the appearance of an elliptical; however the l ight 
profile is best fit by an exponential disc model. Additionally, IHoI 
( l2005h estimate a significant star formation rate (19.4 AlQyr~^) 
based on the [O II] luminosity. Our results indicate a 60.8±27.6% 
intermediate population, a 37.4±I7.7% old population and a sta- 
tistically insignificant young population. Multiple residual narrow 
emission lines are noted for this object. 



REFERENCES 

Bahcall, J. N., Kirhakos, S., Saxe, D. H., & Schneider, D. P. 1997, 

ApJ, 479, 642 
Barthel, P. D. 2006, A&A, 458, 107 
Bell, E. R, et al. 2004, ApJ, 608, 752 

Bennert, N., Canalizo, G., Jungwiert, B., Stockton, A., Schweizer, 

R, Peng, C. Y., & Lacy, M. 2008, ApJ, 677, 846 
Benson, A. J., Bower, R. G., Prenk, C. S., Lacey, C. G., Baugh, 

C. M., & Cole, S. 2003, ApJ, 599, 38 
Bershady, M. A., Andersen, D. R., Verheijen, M. A. W., Westfall, 

K. B., Crawford, S. M., & Swaters, R. A. 2005, ApJS, 156, 311 
Best, P. N., Kauffmann, G., Heckman, T. M., Brinchmann, J., 

Chariot, S., Ivezic, Z., & White, S. D. M. 2005, MNRAS, 362, 

25 

Boroson, T. A., Persson, S. E., & Oke, J. B. 1985, ApJ, 293, 120 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Canalizo, G., & Stockton, A. 2001, ApJ, 555, 719 
Canalizo, G., Stockton, A., Brotherton, M. S., & Lacy, M. 2006, 

New Astronomy Review, 50, 650 
Canalizo, G., Bennert, N., Jungwiert, B., Stockton, A., Schweizer, 

R, Lacy, M., & Peng, C. 2007, ApJ, 669, 801 
CardeUi, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 

245 

Cattaneo, A., Dekel, A., Devriendt, J., Guiderdoni, B., & Blaizot, 
J. 2006, MNRAS, 370, 1651 
Chabrier, G. 2003, PASR 115, 763 

Cid Fernandes, R., Mateus, A., Sodre, L., Stasiiiska, G., & Gomes, 

J. M. 2005, MNRAS, 358, 363 
Condon, J. J., Cotton, W. D., Greisen, E. W., Yin, Q. R, Perley, 

R. A., Taylor, G. B., & Broderick, J. J. 1998, AJ, 115, 1693 
Dunlop, J. S., Taylor, G. L., Hughes, D. H., & Robson, E. I. 1993, 

MNRAS, 264, 455 
Dunlop, J. S., McLure, R. J., Kukula, M. J., Baum, S. A., O'Dea, 

C. P, & Hughes, D. H. 2003, MNRAS, 340, 1095 
Falomo, R., & Ulrich, M.-H. 2000, A&A, 357, 91 



Fukugita, M., Shimasaku, K., & Ichikawa, T. 1995, PASP 107, 
945 

Gebhardt, K., et al. 2000, ApJ, 539, L13 

Govoni, R, Falomo, R., Fasano, G., & Scarpa, R. 2000, A&A, 
353, 507 

Guyon, O., Sanders, D. B., & Stockton, A. N. 2006, New Astron- 
omy Review, 50, 748 

Hamilton, T. S., Casertano, S., & Turnshek, D. A. 2002, ApJ, 576, 
61 

Ho, L. C. 2005, ApJ, 629, 680 

Hooper, E. J., Impey, C. D., & Foltz, C. B. 1997, ApJ, 480, L95 
Hopkins, P. R, Hernquist, L., Cox, T. J., Di Matteo, T, Robertson, 

B. , & Springel, V. 2006, ApJS, 163, 1 

Hopkins, P R, Cox, T. J., Keres, D., & Hernquist, L. 2008, ApJS, 
175, 390 

Hughes, D. H., Kukula, M. J., Dunlop, J. S., & Boroson, T. 2000, 

MNRAS, 316, 204 
Hutchings, J. B., Johnson, I., & Pyke, R. 1988, ApJS, 66, 361 
Hutchings, J. B., & Neff, S. G. 1992, AJ, 104, 1 
Hyvonen, T, Kotilainen, J. K., Falomo, R., Orndahl, E., & Pur- 

simo, T. 2007, A&A, 476, 723 
Jahnke, K., Kuhlbrodt, B., & Wisotzki, L. 2004, MNRAS, 352, 

399 

Jahnke, K., Wisotzki, L., Courbin, R, & Letawe, G. 2007, MN- 
RAS, 378, 23 

Kauffmann, G., et al. 2003, MNRAS, 346, 1055 

Kellermann, K. I., Sramek, R. A., Schmidt, M., Green, R. R, & 
Shaffer, D. B. 1994, AJ, 108, 1163 

Kormendy, J. 2000, arXiv:astro-ph/0007400 

Lacy, M. 2006, Astrophysics Update 2, 195 

Letawe, G., Magain, P., Courbin, R, Jablonka, P., Jahnke, K., 
Meylan, G., & Wisotzki, L. 2007, MNRAS, 378, 83 

Mateus, A., Sodre, L., Cid Pemandes, R., Stasinska, G., 
Schoenell, W., & Gomes, J. M. 2006, MNRAS, 370, 721 

McLure, R. J., Kukula, M. J., Dunlop, J. S., Baum, S. A., O'Dea, 

C. P, & Hughes, D. H. 1999, MNRAS, 308, 377 
Miller, J. S., & Sheinis, A. I. 2003, ApJ, 588, L9 

Nolan, L. A., Dunlop, J. S., Kukula, M. J., Hughes, D. H., Boro- 
son, T, & Jimenez, R. 2001, MNRAS, 323, 308 
Oke, J. B., et al. 1994, SPIE, 2198, 178 

Peletier, R. R, Valentijn, E. A., & Jameson, R. R. 1990, A&A, 233, 
62 

Press, W., Teukolsky, S., Vetterling, W., & Flannery, B. 1992, Nu- 
merical Recipes in C: The Art of Scientific Computing (Cam- 
bridge, MA: Cambridge Univ. Press) 

Schlegel, D. J., Finkbeiner, D. P, & Davis, M. 1998, ApJ, 500, 
525 

Sheinis, A. I. 2002, Ph.D. Thesis, Univ. California 
Tremonti, C. A. 2003, Ph.D. Thesis, Johns Hopkins Univ. 
Vanden Berk, D. E., et al. 2001, AJ, 122, 549 
Vanden Berk, D. E., et al. 2006, AJ, 131, 84 
Wolf, C, et al. 2004, A&A, 421, 913 
Wolf, M. J., & Sheinis, A. I. 2008, AJ, 136, 1587 
Wright, S. C, McHardy, I. M., & Abraham, R. G. 1998, MNRAS, 
295, 799 

Yu, Q., & Tremaine, S. 2002, MNRAS, 335, 965 
Zakamska, N. L., et al. 2003, AJ, 126, 2125 



© 0000 RAS, MNRAS 000, 000-000 



